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Abstract 

The aim of these notes is to show how the magnetic calculus developed in [MP! IIMP11 IIMP21 IMPRI ILMRj 
permits to give a new information on the nature of the coefficients of the expansion of the trace of a function 
of the magnetic Schrodinger operator whose existence was established in |HR2| . 



1 Introduction 

Let us consider the magnetic Schrodinger operator on R defined by 



d 

P A (h) = £(&!>„, - A 3 {x)f + V(x) , (1.1) 

3=1 



where D T :— —id T . and we assume: 

Hypothesis 1.1 

• % £ I c]0, +oo[, with X a bounded set having as accumulation point, 



• A = (Ai, . . . , A d ) with Aj e C c 



• v e c*°°, y > -c. 



It is known that the operator associated with P A (h) on Cg°(R d ) admits a unique selfadjoint extension on 
L 2 (R d ), which can be defined as the Friedrichs extension. We denote by P A (h) this extension. For any function 
g £ Cq°(R), we can define by the abstract functional calculus g[P A (H)). 

We shall also make the following assumption concerning the potential function V: 

Hypothesis 1.2 

S y :=liminfF(x) >infV. 

\x\— >oo 

It is known in this case by Persson's Theorem (see for example |Ag| ) that the spectrum is discrete in ] — oo, Ey [ 
and using the max-min principle one shows easily that the spectrum is non empty for h small enough. In particular, 
for suppg CC] — oo, liminf V[, one can consider Tr g(P A (h)). Our goal is to analyze the expansion of this trace as a 
power series in h and the dependence of the coefficients on the magnetic field, i.e. the two-form B := d(Y]j Ajdxj). 

Of course if we have two vector potentials A and A, such that dA = dA = B, we know that there exists <j> £ C°° (M. d ) 
such that : A = A + dcj>, and the conjugation by the multiplication operator by exp 40 gives a unitary equivalence 

between P A (h) and P A {K). Hence Tr g(P A {Kj) and its expansion should depend only on B. We would like to 
investigate how it depends effectively on B. 

Our main theorem is the following. 

Theorem 1.3 _ 

Under the previous assumptions on A and V and with H = P A {K), there exists a sequence of distributions T B £ 
Z>'(R), (j £ N) such that for any g £ C^°(R) with supp .<? c] — oo, Ey[ and for any N £ N, there exist Cn and Hn, 
such that: 



{2nh) d Tr g{H) - Y. hJT f(9) 

0<j<N 



< c N h N+1 , \fh e]o, h N ] n I. (1.2) 



More precisely there exists kj £ N and universal polynomials Pi(u a , vp,j t k) depending on a finite number of variables, 
indexed by a £ N 2d and f3 £ N d , such that the distributions: 

T f(9)= E [g {i) (F(x,0)Pdd^F(x,0,de x B jk (x))dxdt;, (1.3) 



0<«<fe 



w/iere F(a;,0 = £ 2 + V(x), safe/y ([T2]) . Finally, Tf = /or j odrf. 

This theorem was obtained under stronger assumptions in [HRlj . but the main difference with the statement 
above was that the expression of T^(g) was given in terms of a vector potential A such that dA = B. Tricky 
calculations permitted after to recover a gauge invariant expression for the three first terms : 

T B (g):= [dxd£g(F(x,0), T 1 B (g):=0, T?(g) := -1 / dxdtg"(F(x,£)) [(AV)(x) + \\B(x)f] . (1.4) 



The approach of [HRlj did not permit to recover the same kind of result for any term of the expansion. On the 
contrary, we will show that, when it can be applied the magnetic calculus permits to give naturally this expression. 
To state the results at the intersection of the domains of validity of the two calculi is actually unnecessary. Following 
essentially arguments presented in IHMRj in the case without magnetic potential, we will show how we can use the 
Agmon exponential decay estimates in order to modify the behaviour of V and A at infinity, without changing the 

asymptotic behavior of Tr g(P A (h)), in order to enter simultaneously in Helffer- Robert's class and in the magnetic 
pseudodifferential calculus of [IMP 11 HMP21 iMPllMPR] , 

The paper is organized as follows. In Section [2j we review the now standard ?i-pseudodifferential Weyl calculus 
and the attached functional calculus. Section [3] is devoted to the presentation of the gauge-invariant magnetic 
calculus [IMP 11 IIMP2J [MPllMPRj . The last section will give the proof of the main theorem. 

We shall constantly use the notations X = M. d , 5 := X x X' , with X' the dual of X. The points of H 
will be denoted as X = [x, £). Recall that H has a canonical symplectic form a((x,£,),(y,r])) := £(y) — rj(x) = 

S (ZjVj ~ Vj x j)- We denote by C^(y) the space of C°°-functions on the vector space y having at most 
i<j<d 

polynomial growth at infinity together with all their derivatives. We denote by Cf^, u (3^) the subspace of functions 
having a unique polynomial upper bound for all their derivatives. 

Acknowledgements We ackowledge the support of the CNRS Franco- Romanian European Associated Labora- 
tory Math-Mode and of the University Paris Sud as well as the ANCS Partnership Contract No. 62-056/2008. The 
second author thanks the University Paris Sud for its kind hospitality during part of the elaboration of this paper. 

2 The /i-pseudodifferential Weyl calculus and trace formulas 

To any Schwartz test function <j) G 5(S) we associate a bounded linear operator Op h ((f>) on the Hilbert space 
H:=L 2 (X): 

[Op h (4>)u] (a;) := (2nh)- d J^J^dyd V erM-y) ( j ) (^-,v) «(»), Vu G 5(H), 
for some constant h > 0. It is easy to prove that Dp h (<fi) G M(H) and 



\\Op h (4>)h(.H) < / dx\\[$-<f>](x,-)\\oo 
Jx 

where $~ denotes the inverse Fourier transform with respect to the second variable. Moreover it is not hard 
to prove, by using Schur's Lemma, that Op h : 5(S) — ► M(H) extends to an isomorphism of topological vector 
spaces Dp^ : S'(E) — > M(S(X);S'(X)). We can transport the operator multiplication from B(H) back to a non- 
commutative product on 5(S) 

DpMOPkW =■■ OpnMitf), Vfa, V) G [5(S)] 2 . 
Explicitely we have 

($M0 (A) := (nh)- 2d f dY J dZ exp [- (2i/K) a(Y, Z)\ cf>(X - Y)^{X - Z). (2.1) 

One can prove that 

Dp fe [5(3)] =M(S(X);S(X))=M(S'(X);S'(X)) 

so that one can consider products of the form Dp h (^)Op h ((t>) and Dp h ((/))£)p h (®) for ($,</>) £ <5'(S) x 5(H) and 
define the Moyal algebra 

m(E) := {$ e S'(Z) | $(W G S(E),<M h $ g 5(H), V^ G 5(H)} , 

that is obviously an algebra for the tj^-multiplication and even a ^-algebra for the anti-involution given by complex 
conjugation of distributions. 

The limit h —$■ 0, that has a rather singular behaviour, should correspond in some sense to the classical algebra 
of observables that is a commutative algebra. A precise meaning of this limit can be given in the context of strict 
deformation quantization (see |Laj ). but is out of our aims in this paper. On the contrary, the asymptotics K — > 0, 
known as the semi-classical asymptotics, and considered in the frame of asymptotic series in H is an important 
problem and we shall concentrate on some of its aspects when magnetic fields are present. 



The Moyal algebra contains many interesting subalgebras, among them the usual Hormander symbols 
S?(X) := IfeC 00 ^) I sup <£>-™+M \(dZd?F)(x,£)\ < oo I . 

1 (x.f)eH " J 

A symbol F G Si*(X) of strictly positive order to > satisfying 

3C > 0, 3R > such that C < £ > m < |F(a;,f)|, V(a;,£) G 5 with |x| > i? 

is called elliptic and has the property that it exists a positive constant a > such that for any j £C\lU (oo, —a) 
the distribution 3 + F has an inverse (3 + F)~ with respect to the (U-product and this inverse belongs to the class 
Si m (X). In other words, the operator Dp h (F) has a self-adjoint extension with a resolvent that has a symbol of 
Hormander class S^ m (X). An important problem is how to relate this pseudodifferential calculus (defined by the 
Moyal product jj^) with the usual functional calculus for self-adjoint operators, when these operators are of the 
form Op h (F) with F G S?(X) elliptic (to > 0). 

In dealing with semi-classical problems, the parameter h is no longer constant and it is important to consider 
asymptotic series in H. An essential fact is that the Moyal product U has a 'suitable' behaviour with respect to 
such asymptotic series. More precisely, let us consider the space S[ s ' m (X) of ^-symbols of the form F : X x 5 — > C 
such that 

F(h)eC°°{~), \/hei, sup h- s < £ >- m +H | (d^F) {h, f, x) | < 00. 

IxH 

We refer to |RobllHe2] for a systematic discussion and for more general ^-symbols. In the frame of these symbol 
classes, the Moyal product U has the following property (see |RobMHe2] ): 

V(F,G) G s[ Sl ' mi) (X) x s[ S2 ' m2) (X), FUG G ^«+'"'»i+ TO »)(^ and 

FUG -FGe s^+^+^+^-^iX). 
We shall usually work with elements F G S^ (A") having an asymptotic expansion of the form: 

F(h,^x)^h s J2^ k Fk(x,0, withF k ESr- k (X). (2.2) 

keN 

The symbol h s Fo G 5'} s ' rra (X) is then called i/ie principal symbol of F. 

What is important here is to find a class of functions (actually essentially C§°) such that g(F) is a nice pseudo- 
differential operator with simple rules of computation for the principal symbol. We are starting from the general 
Dynkin-Hclffer-Sjostrand formula (see |DiSj| ) 

g{P) = -it- 1 lim // |r(A,/*) (A + i/x- P) _1 dAd/i, (2.3) 

which is true for any selfadjoint operator P and any g in C(j° (M) . 

Here the function (A, /x) 1— >• g(A, fi) is a compactly supported, almost analytic extension of g to C. This means that 

g = g on R and that for any N G N there exists a constant Cat such that | gf(A, /x)| < CatI/xJ^ . 

The main result due to Helffer-Robert (see also [DiSj] and references therein) is that, for P = D\) h (F) a self-adjoint 

operator with a /i-symbol F real and semibounded from below and having an asymptotic expansion as above (J2.2I) 

with s = and g in Cq°(M), the operator g(P) is a ft-pscudodiffcrcntial operator of the form Dp h {gh(F)), whose 

Weyl symbol g~fi(F)(h,£,x) admits a formal asymptotic expansion in h 



gH(F)(h,Z,x)~y]h k g k (F)(x,0, (2.4) 

with 



fe>0 



g (F) =g(F ), 

9l (F) =F 1 -g'(F ), (2.5) 



^2fe-l 



9k(F) =EZid k , l g^(F ),Vk>2, 



^1 



where the dk,i are universal polynomial functions of the symbols d^d^Fg , with |a| + \/3\ + £ < k . For h = 2, one 
has 

*,! = ^2, ^2,2 = P2,2 + (1/2)*?, d 2 , 3 = P2,3 (2.6) 



with 

t ^ly/ j^Q & f q d2p Q d2p Q \ f27) 

P2,2{X,S) 8 2^y d d^dXkdtj dXjdXkdZjdtJ { -'> 

V2t(x, = y a E ( 2 C^o ■ d *i F ° ■ d Z* F ° ~ <-*^o • d^Fo ■ d ik F - dl ik F ■ d Xj F Q • d Xk F ) . (2.8) 

j,* 

The main point in the proof is that one can construct for Qz /Oa parametrix (= approximate inverse) for (P — z) 
with a nice control as 3z tends to 0. The constants controlling the estimates on the symbols are exploding as 
3z — » but the choice of the almost analytic extension of / absorbs any negative power of |5z|. 
As a consequence, one gets that for H small enough, if for some interval / and some eo > 0, 

Fq 1 ^ + [— eo, eo]) is compact, (2.9) 

then the spectrum of Dp h (F ) is discrete in /. In particular, one gets that, if Fq(x,£) — > +oo as \x\ + |£| — > +oo, 
then the spectrum of Op h (F ) is discrete (Op h (F ) has compact resolvent). In fact one gets more precisely the 
following theorem (due to Helffer- Robert). 

Theorem 2.1 . 

Let P = Qp h (F) be a self-adjoint operator with a h-symbol F real and semibounded from below, having an asymptotic 
expansion of the form \2. 2)) and also satisfying V2. 9]) with I — [Ei,E2], then, for any g in Cq°(]Ei, E2D, we have 
the following expansion in powers of h : 



Tr [g{Dp h (F))] ~ pirh)-*^ T M , (2-10) 

j>0 



where g *— > Tj(g) are distributions in D'(]Ei,E2[). 
In particular we have, when F\ = F2 = 0, 



Mg) =JJg(F (x,O)dxdi, 

T iG?) = ' (2.11) 

Ug) =-kH9"{F {x,0)Y, j , k Sfeafe - l&- k &k) *"*- 



This theorem is obtained by integration of the symbol of g(Op h (F)^j given in (J2.4I) . because we have the needed 
regularity so that the trace of a trace-class operator Dp h (F) is given by the integral of the symbol F over 5. 
According to the definition of the Weyl quantization, the distribution kernel is given by the oscillatory integral: 

K(x, y] h) = (2nh)- d J^exp^(x- y )-^j F (h, £, ^±l) d£ , (2.12) 

and the trace of Dp h (F) is the integral over X of the restriction to the diagonal of the distribution kernel: K(x, x) — 

(27rh)- d J x F(h^,x)di. 

Of course, one could think of using the theorem with g the characteristic function of an interval, in order to get 
for example, the behavior of the counting function attached to this interval. This is of course not directly possible 
and this will be only obtained through Tauberian theorems ( [Holj . |Ho4j and [Iv]) and at the price of additional 
errors. Let us however remark that, if the function g is not regular, then the length of the expansion depends on 
the regularity of g. So it will not be surprising that, by looking at the Riesz means: g s ,E(t) '■— max{0, (E — t)} s 
(for some s > and E £ (E±, £2)), we shall get a better expansion when s is large. 

Under some assumptions on A and V, including the condition divA = 0, one can show that Pa(K) is an H- 
pseudodifferential operator whose total Weyl-symbol is F(x, £) = (£ — A) 2 +V in some class of [HRlj ■ More prcsely, 
we have to assume that for any a £ N d , we have 

\d%V{x)\<C a {V(x)+C + l), (2.13) 

and, for j = 1 . . . , d, the following non gauge covariant conditions: 

|^A,(aO| <C a (V(x)+C+ 1)3. (2.14) 



3 The /i-magnetic quantization 

3.1 Results for fixed h 

We consider a magnetic field described by a bounded smooth closed 2-form B on X = M. d and the associated 
modified symplectic form on 2 

a (x,o((Vi V), (C, z)) := r)(z) - ((y) + B x (y, z), 

that may be used to define the classical Hamiltonian system in the given magnetic field. For the quantum description 
we shall also consider the Hilbert space T-L = L 2 (X), we shall choose a smooth vector potential A, i.e. a 1-form 
satisfying the equality B = dA and we shall define the following gauge covariant representation, for all <j> G 5(H) 
and all u G 5(H), 

'Op£(cf))u\ (x) : (2Trhy d [ J dydrie^-^e-^^.y^cfrf^-^,^ u{y) , 

where J, ,A denotes the integration of the 1-form A along the oriented segment [x,y]. This gauge covariant 
'magnetic quantization' allows to define a 'magnetic' Moyal product |}j? 

(#f t/>) (X) := (nh)- 2d J s dY f s dZ exp [- (2i/h) o* (Y, Z)] <j>{X - Y)yj(X - Z) 

= (nh)- 2d J s dY J s dZ e -{^m<T(Y,z) e -(i/n)6 B (x, v ,z^ x _ Y)$(X - Z) (3.1) 

= (irh)- 2d jldY jldZ e -(2im<r(x-Y,X-z) e -(i/h)<> B (*,v,z) ( f ) (Y)ip(Z). 



6*{x,y,z):= I B , 9 B (x,y,z):= / B 

J <x— y— z,x-\-y— z,x— y-\-z~> J <x — y+z.y — z+x,z— x-\-y> 

Here < x — y — z,x + y — z,x — y + z > denotes the triangle defined by the three points x — y — z,x + y — z, and 
x — y + z, with the usual trigonometric orientation and the integrals of B denote the integrals of the two form on 
the given oriented triangle. Associated with this product we can define a 'magnetic' Moyal algebra for the magnetic 
field B: 

M B (E) := {$ G 5'(H) | nU G 5(H), n B $ G 5(H), V0 G 5(H)} , 

This 'magnetic' Moyal calculus preserves a large number of the nice features of the usual Moyal calculus and we 
shall recall some of them that are useful for our analysis of semiclassical trace formulas. 

Proposition 3.1 (Propositions 3.5 and 3.10 in \MFp 

For any magnetic field B with components of class C^(X), one can find a vector potential A with components 
also of class C^(X) and then the application Op^ : 5(H) — > M(L 2 (X)) extends to an isomorphism of vector spaces 
Qp£ : 5'(H) -> B(S(X);S'(X)). The above isomorphism has a restriction Dp£ : L 2 (H) -> M 2 (L 2 (X)) that is 
unitary (here 1S>2(L 2 (X)) is the algebra of Hilbert- Schmidt operators on L 2 {X)). 

Proposition 3.2 (Proposition 4.23 in JMP\j and Lemma 1.2 in flMPlf ) 

For any magnetic field B with components of class C^ (X) , we have the following inclusions: 

C~ liU (H) C OT B (H); S?(X) C M B (E). 

Proposition 3.3 (Theorem 2.11 in \LMR^ ) 

For any magnetic field B with components of class BC°°(X) the 'magnetic' Moyal product defines a continuous 

map 

S™ 1 (X) x S™ 2 (X) 3 (F, G) *-+ F)jf G e S{ ni+rn2 (X), 

and for any N G N there is a canonical expansion 

iV-l 

Fjjf G =J2 H J + r n, with H,j G S? 1+m2 - j {X), R N G S™ 1+m2 - N (X) 

3=0 

in which Hn = F ■ G. 



Proposition 3.4 (Theorem 3.1 in JIMPlf ) 

For any magnetic field B with components of class BC°°(X) we have that for any associated vector potential A, 

£>*£[£?(*)] cl(£ 2 W) 

and there exist two positive constants c, p depending only on the dimension d of X such that 
DP ^ F) , ((i(w ^ su P( su P( SU P <Z> H \(d a x d?F)(x,0\)). 

HL 2 (X)) | |< p | a |< p ( X ,£)£H 

Proposition 3.5 (Theorems 4.1 and 4.3 in UMPlf and Proposition 6.31 in UMP2f ) 
Suppose the magnetic field B is of class BC°°(X); then 

• if F £ Si(X) is a real function, 0\) h (F) is a bounded self-adjoint operator on L 2 (X) for any vector potential 
A of B; then the resolvent of Dp h (F) has a 'magnetic' symbol of class Sf(X); 

• if F £ S™(X) is a real elliptic symbol with m > 0, then Op h (F) has a self-adjoint extension in L 2 (X) for 
any vector potential A of B and the resolvent has a 'magnetic' symbol of class Si m (X); if we choose A with 
components of class C^(X) then D\) h (F) is essentially self-adjoint on S(X) and its self-adjoint extension 
has as domain the 'magnetic' Sobolev space: 



U; 



XX) := [u £ L 2 (X) | Dp£(p m )u £ L 2 (X), where p m (x,£) :=< £ >'"} ; 



• if F £ S" l (X), with m£i, satisfies 3lF(x,£) > C|£| m for |£| > R, with some strictly positive constants C 
and R, then for any r el there exist two positive constants Co and C\ such that for any u £ W^ [X) we 
have 

5ft < u,Op A (F)u > L 2 (X > C \\u\\ h a /2(x) - dWuWu^x) . 

Proposition 3.6 (Proposition 6.33 in \IMP2j ) 

Suppose the magnetic field B has components of class BC°°(X) and $ £ Cq°(K), then, for any real F £ SY l (X) 
with m > 0, F elliptic for m > 0, and for any A such that dA = B, the operator $>(Dp h (F)) , defined by the 
functional calculus for self-adjoint operators has a 'magnetic' symbol of class S^ m (X). 

3.2 Semiclassical results 

Let us consider now the dependence on H £ I. We shall come back to (|3.ip and analyze the ^-dependence of this 
product: 

(#jfy) (X) = (irh)- 2d I dY J dZe- ( - 2i /^ a{ - Y ' z ^e- { - i ' n ^ B{ - x ' v '^4>{X - Y)ij){X - Z) 

with 



3 B 



(x,y,z):= f B = A^y jZk f ds f dtB jk (x + (2s - l)y + (2i - l)z). 

J <x-y — z,x+y — z,x — y+z> -j, JO JO 



Let us make the change of variables (y, z) t-t (Hy, Hz) in order to obtain 

mU){X) (3.2) 

= Tr- 2d I o\yI dZ e -iMv,^ c -(<>9** (».iM)0(a; - %, £ - rf)ip(x -hz,£-Q, 

with 

e%{x,y,z)=A y Vy j z k I ds [ dtB jk (x + (2s - 1)% + (2t - l)Hz). 

We notice that we can now obtain an asymptotic expansion of the ^-product with respect to h by using the 
Taylor formulas for (/>, tj;, like in the non-magnetic case, and for the exponential e~ i -- lh ^ en ^^'^ and also for B in 
the expression of Q%(x, y, z): 

4>{x-hy,t;-ri)= ]T tK £ ^(^(^ _ „) +^ Nj 

0<v<N ' \a\=v 



i,(x-nz^-o= E ^E^(^)(^-o+5K 



4>,N, 



0<n<N 



l/3|=M 



r (*)tf (.,„.) = ^ L^ [#(*,,,, 3)]' + *^, 



0<p<W 



(a;,j/,z)= E TT 



0<A<N 



E E ^ z * ( a7 ^ fe ) (*) I ^ / ds [ dt \ s y + iz ) 

|7| = A \jk J T - 7 - 1 - 7 - 1 



■*B.N 



^ A! 

0<A<Af 

where for any TV > 1 in N we have 



EE r ^^(E^( 97 ^)( 3 



7| = A (5<7 



jk 



+ *B,N, 



yCt pi ^ __ 

<R^ N (h,x,y,£„rj) = {-K) N+1 E ~ ( d x 4>) (x - uhy, £ - T])du = E V a ^4>,N,a(^x,y,^,r]), 

\a\=N+l ' ° |a|=JV+l 



|/9|=JV+1 



|/3|=JV+1 



iJV+1 



Jo Jo 



ft + XN+i(ft,x,y,z) 



E 

0<A<iV 



A! 



EE^y^[E^(^)( 3 



■7 ,,<5 r 7-<5 
|7|=A 5<7 y jk 



*B,N 



N+l 



?B,N(ft, x,y, z) = ^ JV+1 E (tO 7 ^ s /^ / duy^^yjZk (d 1 Bjk) {x + ftu(sy + tz))(sy + tz) 1 
| 7 |=JV+i ^ ^- lJo jk 

= ^ +1 X;^fcE (7!)- 1 ^, fc (^^2/,^)E T I ^*- 

jfc | 7 |=jV+l (5<7 

Let us observe first that the powers y a and z@ appearing in the Taylor expansions of <j> and i/j are dealt by integration 
by parts, using the 'oscillatory' exponential e ~ 2la '( Y - z ) an d transformed in ^-derivatives of ip and respectively in 
^-derivatives of 4>. A similar procedure may be used to get rid of the factors y z 1 appearing in the expansions of 
6f*(x, y, z) p ; in fact we have 

j," e -«"( y .*) = H^H^e" 2 *^; zee- 2 *™ = (ifflWtfe- 2 *™. 
Thus, denoting by 



ii \ Ik 
|7|=A ■> 



one has 



*?(*):= E E(- 1 ) AH5lT ?( 97 ^)w 5 c% 5 r^ 

5<7 
(^ V)(X) = r.-- -'' / r/V/ ,'Z ( - 2 ' Vrn - Z) X 

(-ift) p 



)//,- - 



(3.3) 



^ p! 

0<p<N ' 



ft X fl^ X 



0<A<JV ' v 



+ 5Hs.7V> X 



E 
E 






A 



E ^(WW^t^-^+iR 



0,JV 



0<1/<JV \a\ = v 



(-ft) 1 " v- /" ! 



0</j<W 



l/3|=M 



(2tt)- m £ 



-ihy 



0<p<N 



A / „• \ A 



(j)^+P(_l)P(^)^+M 



^ 4A! \2 

0<A<W 



/! ill 



£?(*) 



(3.4) 



E w -^,7 - E 2gWw)(..o(flrW(-.o 



0<^<JV 



|a|— iv 



+^" 2d / dW dZ e -**V,z) [ tB)JV ( S) a . j yj Z )3 N {H, x, y, z) + d\ BiN (h, x, y, z)\ <f>(x - %, £ - rj)1>{x -hz,£- C) + 



i \ 
2 



-2d 



N+l 



dY dZ e - 2i °™ { £ 



(-iny 



0<p<N 



E 

0<A<7V 



h x /*\ A 



4A! V2 



£*(*) 



/ I -. I 7V7 I 1 



-(-1) 



\a\=N+l 



N+l I l 



2(JV+1) 



E E (^^,^a)(M,y,^£,?/)(d"%/,,A^)(ft,2:,2/,^£,C) 



a|=AT+l |/3|=JV+1 



Let us discuss the remainders in the above expansion. First let us consider the factor 3Jv(A x, y, z) 

h x /r* 



3 N (h,x,y,z)= Yl — r~ E C p x b,n%x,v,z) 

l<p<N "' 0<K<p-l 



p — K 



0<\<N 



^k& «<•> 



and taking into account the differential operators contained in lf(x) (see (|3.3p ) let us compute for p > 1 and 
< K < p — 1 one of the terms appearing in the above sum (with < Xj < N for 1 < j < k) 

drf dZ e- 2l ^ Y - z \ B , N {h, x, y, zy- K ll (x) ■ ■ ■ %f H (x)<f>(x - hy, £ - r,)^{x -hz,£- Q 



= ftt^X"—) / dY dZ e - 2l ° (Y ' Z *> 



E^ Zfc E W)- 1 F% j ,k(?i,x,V,z)'52 T 8 v : 



1 ,,<5 r 7-<5 



jk |7|=JV+1 



<5<7 



p—K 



xT£ (or) • • • Zf H (x)4>(x - %, £ - r?)V(a; - &*, £ - 



AT+l 



h (N+l)(p-K) 



£ / a «/ s ^- 2, '™0«. W )* 



l"l=l^l=p-« ' 
x ^(-l^+MSlT^-^Cz) . • .£?»(^)(:c - fiy,£- r,)(d^)(x -hz,£- C) (3.5) 

<5<7 

where G^o is a product of p — k functions of class BC°°(X 3 ) uniformly for h e (0, Ho], depending only on the 



I 5 Ac q™i-|/3| 



V*.l, a Q™2-|a|, 



derivatives of order N + 1 of the magnetic field B. Moreover 9? cf> <E S x 1 ' \X) and dfrj) £ S 1 2 (X) uniformly 

in h G (0, Hq] so that the integral (|3.5[) defines an element in S{ 
Let us consider the second contribution: 



(JV+l,mi+m 2 -2-(JV+l)), 



(X) for any p — k > 1. 
dW dZ e" 2 '" T(y ' z) 5n B ,7v(?i, aj, y, z)<j>{x - %, £ - ij)tp(x -hz 7 £_- () (3.6) 

= {-ih) N+1 (jdui... f UN e- iUN +^ e S^y^du N+1 ^j x 

N+l 



x / dW dZe- 2 ^ y ^[^(x,2/,z))] + <f>(x-hu,£-ri)1>(x-hz,S-C) 



so that a similar procedure with the one used for (|3.5[) proves that this integral defines a function of class 

o(JV+l,mi+m 3 -2(iV+l))/«,v 

For any </> 6 S m {X) the rest iSH^jv.a is a function of class 5} ,m ~ ( (<^0) Va £ N d so that it is easy to 

notice that the last contribution to the rest: 

dYf dZ e - 2i ^ Y ' z U J2 



-2d 



{-my 



0<p<N 



0<\<N 



H Hi U ^ 



) E (^,iv, Q (a c Q v)(x-fe,e-C) + (-i) Ar+1 ^0,^a(9 I »(s~%,e-»/)) 

' ~. I — AT I 1 



2 



q|=AT+1 



JV+1 / * 



2(JV+1) 



^ 1)JV+1 (i) ^ ^ (a?^,a)(fi,»,»^^,»?)(^'^,^/s)(ft,a!, y ,«,e,o 

^ ' |a|=iV+l |/3|=iV+l 

defines an element of sf+^+^-^+D) {x y 

Let us now concentrate on the main terms in the expansion (|3.4p : 

p 



-Id 



dY dZ e^< Y ^ \ J2 



{-my 



0<p<N 



pi 



E 

0<A<Af 



4A! V2 



£?(*) 



{-ny 



0<u<N 



E ^rESwcr^K^-^) 



(-^ 



E ^Es(Hm/(^)(^-o 



= 7r- 2d / dW dZ e - 2 -(^- z ) J ^ 



H*)' 



0<H<N |/3|=M 

Ai+-+A p 



/3! 



4Pp\ 

0<Aj<p 
l<j<p 



£ (t 



Ai!---A p 



■^(^•••Tffx) x 



i/+/j 



0<i/<JV 
0<)J,<N 



E (-ir(?) E ^(^)(*.*-*)(^)(*.e-o 



|a|— i/ 



E ^ k E E 



a!/3 



(_!)P+i2jP 



E 



0<fc<JV 



•\ d+b + Ai + .-.+Ap -^ 



l +l 2 +l=k p+\i+...+\„=l 

a<h<k o<p<i 

Q<l 2 <k 0<Aj<i 

0<l<k l<j<P 



AP P \ ^ \ 2 

I0M 



o 



\R\ 



(3.7) 



E ••• E {-lf ll+ - +M Tl^--T^{d^B llkl ){x)---{d^B Jpkp ){x)x 



l7l|=Ai 
(5i<7i 



|7pI=A p 
<5 p <7p 



x (% x •••% p ^ +7l ^ 1+ '" +7p ^9^) (s,0 (d( il -d iip ^ +Sl+ - +s "^i>) {x,0 + R N {h;<t>^){x,0, 
where i? w £ ^+1^+™^+!)) ( ^ 

Using these expansions one gets the following statement. 

Proposition 3.7 

// i/ie magnetic field B has components of class BC°°{X) and <p £ S T " 1 {X), ip £ S™ 2 {X) then there exists a 
sequence {c^(0, -0)}feGN such that cf{cp,ip) £ c™i+m 2 - (^ an ^ j or an y TV > 1 m N, we have: 



--N 



N-l \ 

fe=0 / 



+m 2 -iV 



(*) 



uniformly for liGl. For /c = 0, we /save Cq {(f>,ip) = (j> ■ tp. Moreover, for any k £ N, £/ie function C?{(j>, if); X) only 
depends on values of the magnetic field B and its derivatives of order < j — 1 evaluated at x. 
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We shall use the notation 

N-i 
7f(«^):=#fV- <f>-i>=Y,h k c%(ct>^)+h N i%(<f>,i;), (3.8) 

fc=i 

witht^,^) G s{°^ 1+m2 - N) (X). 

Starting from the above result we shall work in the sequel with functions of class 5j s,m (X), thus admiting 
asymptotic expansions in h G I of the form given in the above proposition. 

4 The magnetic Schrodinger operator 

4.1 Preliminaries 

Let us notice that the magnetic Schrodinger operator defined in (|1.1[) satisfies 

P A (h)=Qp£(e + V),. (4.1) 

We shall suppose that the magnetic field has components of class BC°°(X), that the vector potential has been 
chosen of class C^(X) and that V is a real BC°°(X) function. Hence F := £ 2 + V G Sf(X) and all the results 
in Section [3] can be applied in order to conclude that (taking also into account that V is a bounded self-adjoint 
perturbation of P^(H) '■= £>p h (£ 2 ) and the Neumann series expansion of the resolvent) : 

Proposition 4.1 

1. P (H) defined in J4-l\ ) is essentially self-adjoint on S(X) and its self-adjoint extension, denoted H , has the 
domain TL^iX). 

2. The resolvent (H — i)~ l is well defined for 3 G C \ {y G R | y > a} for some a G K and has the form 
(H - 3)- 1 = £>p£(rf ) with rf G si°'~ 2 \x). 

3. There exists do > o depending on B and V such that a ess (H) C [do, 00). 

4.2 The resolvent 

We shall concentrate now on the asymptotic expansion of the symbol rf G S{ ' (X) with respect to ft G I. In 
fact, as the case witout magnetic field is well-known, we shall only be interested in the 'magnetic' contribution to 
the terms of the ft-asymptotic expansion, specifically in putting them in a manifestly gauge invariant form. For 
this purpose we shall use our parametrix type construction of |MPRj in order to express rf in terms of (F — 3) -1 . 
In fact we know that (F — $)^r? = 1. In order to shorten our notations we shall denote by p := F — 3. Let us 
compute 

JV-l 

pMp; 1 = 1 + 7s (p^p; 1 ) = 1 + E ^cfip^p^ + ^i^p^p; 1 ) 



fe=i 

/, Vizi's i c ui v"-; aiiu '■N\P i 'P^ ' v ' ? 

rf=p- 1 -rfttf 7 |(p J ,p- 1 ) 



with cfip^pT 1 ) G V(A0 and t&fo, p" 1 ) G s{°^ N \x). Thus 



0<i<M-l 



E (-i)V»f 

0<i<M— 1 



7Vi-l 1'ftJ 

E^fC^O + ^^'O 

. fe=i 



+(-l) M « 



tv 2 -i 



E^f^.O + ^fc-O 



fe=i 



C« 
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Let us first study the corrections c^(p bl p 1 ). We use the general formulae p.8[) . (|3.7[) and obtain: 



cZ( Ps , P ;i) = (2n)-™ Y, E 



(-1)p+ 1 HP 



h+h+l=k p+Ai + ...+A p =i 
0</i<fe 0<p<l 

0<l 2 <k 0<Xj<l 

0<l<k l<j<P 



|a|=i 

I0H 



1 



a\/3 



1/9! 



£ ••• E (-i) |5ll+ - +l ^ | T7 i 1 ---r7;(^i?, lfel )(x)...(^i? Jpfcp )( a; )x 



7i|=Ai 
<5i<7i 



|7pI = A p 
<5p<7p 



An important observation is that all the derivatives of the symbols p and p _1 have a very special dependence on 
3. More precisely we have: 

• (3?P,)(s,0 = (^)M = (9?V)(*)5 

• (%Pi)(x,0 = (%^)(^0 = 2^; (%%P,)(&,0 = 2<5, fc ; (afp,)(s,0 = 0, V|a| > 3 ; 

• (d^p } )(x,0 = 0, for |a| > 1, |/3| > 1 . 
Lemma 4.2 for any multiindices a and ft we have that: 

(^p,)(*,o= E ifcfoOpr 1- *^) 

0<fe<|a| + |/3| 

where q k (x, £) are polynomials of degree at most k in £ with coefficients functions of x depending only on the first 
\a\ derivatives ofV(x). 

Proof. 

In fact we have: 

(d Xj p- l )(x,0 = -p- 2 (x,0(9 Xj F)(x,0 = -p- 2 (x,0{d Xj V)(x); 

(dt i p; 1 )(x,t) = -p?(x,t)(de J F)(x,t) = -2Z J PT*(*,£)- 

Thus the statement of the Lemma is true for \a\ + |/3| = 1 and we shall proceed by induction on \a\ + \(3\. Suppose 

the statement has been proved for lal + |/3| = N and let us compute the next derivatives. 



(d X] d:d^ Pi )(x,o = d X] 



E 

0<k<N 



IfcP, 



-1-fe 



(*,0 



= E (d Xj q k )(x,0 P ; 1 - k (x,0-(i + k) Yl ^^Oid^^Op; 1 -^^,!;) 

0<k<N 0<k<N 

= E [(d Xj c i k)(x,0-kqk-i(x,0{d Xi V)(x)}p; 1 - k (x,0 - (1 + N)q N (x^)(d X3 V)(x)p; 1 ^ N+1 \x,0; 



0<k<N 



(d^d^ Pi )(x,0 = d^ 



„-l-As 



foO 



E i*p, 1_ 

Cfc<JV 

E (o ij c ik )(x,op: l - k (x,0'^+k) y ^M&^&zK'-^&o 

0<k<N 0<k<N 

E [{dt i q k )(x,Q-2kS J q k - 1 {x,Q]p; 1 - k (x,Z) - 2(1 + N)^ N {x,i)p- l -^ +1 \x, 0- 

N + 1 with coefficients having obviously the same structure. 



0<fc<AT 

Thus the formula is valid also for la 
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Remark 4.3 Some simple computation proves that: 

cf(p J> p- 1 ) = 0, 



'(P,.?,- 1 ) = k"" a (s,0(AV)(s) ~ k" 3 (^0 l(W)(x)| 



+ \p- 2 { Xl \B(x)\ 2 - 2 P ; 3 (x,t) ]T B jk (x)B jm (x){ k U + ^- a (i,OE( fl i B i*) (*)&" 

j fern jfc 

Developing successively each (if -product in (J4.2D and using Proposition 13. 71 one gets the following statement. 

Proposition 4.4 

The 'magnetic ' symbol rf* admits for any N £ N an asymptotic expansion in h of the form 

rf(X)= P ; 1 + J2 Vxf(y ) X) + h N ^;h,X) 

1<J<N-1 

where the terms xf (3; X) only depend on the magnetic field B and its derivatives up to order j — 1 evaluated at X , 
and the rest t^(y,H,X) only depends on the magnetic field B (in a non-local way). Moreover xf '(3) G S l J {X) and 
5(3) € S[°'- N \X). 
Considering (|4.2j) we obtain the following expansion of rf in powers of h 



„B 



P? + Y, W E^ 1 )" E P:HM 1 {P^ 1 )fn--4MM^ 1 ) (4-3) 



Kj Kk<j , . . , 

l<A,<j 

i<i<k 



and developing further each fi-dependcnt jj^-product, one has 



„B 



^+E E(-D fc E E ft ,+w+ - +w C.*(pr 1 .^(p..O.-.^CP„0) 



l<A;<j l<l<k 

l<l<k 

=p; 1 + ^ n E E(-D fc E E £L...»M 1 >^p; 1 )>---^>p; 1 )) 

Kn Kj<n l<k<j \ , , \ ■ , 

1<A,<j 0<« 

i<;<fc i<;<fc 

where 

£ m*.-mi (/' 5i, • • • , fffe) := c* (c* _, (. . . c* (/, 31), 52), ■ ■ ■ , 9k) 
Putting together all these formulae we conclude that: 

Proposition 4.5 

Each term tf (}) for j > 1 is a finite sum of the form 

«?&)= E £(*>Oi>r 2 - p (s,0 

0<p<j 

where fS(x, £) are polynomials in £ of degree at most p whose coefficients are C°° functions of x depending only 
a finite number of partial derivatives of V and B at the given point x. 
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Remark 4.6 Some tedious computation gives: 
^f(3) = -P7 1 cf(p 3 ,pr 1 )-0 ; 

t?(a) = -p 3 _1 cf OvO +p 3 _1 cf {p t ,p- x )ti (p^p; 1 ) - cf (p i _1 ,cf (p.,,?: 1 )) 

= -P t _1 cf0 8 ,P 3 _1 ) 

= -^(x.OCAV)^) + K 4 (^0 |(V^)(x)| 2 

+2 Pj - 4 ( a! ,o E* m (i - ¥im) (d xt d Xm v)(xM m 

-\pf{x, \B(x)\ 2 + 2p- 4 (x, T, jkm B jk {x)B jm {x)^ m - \pf{x, £ jfe (^%0 (*)& 



i3 



tffo) = -P; 1 ci{p il p; 1 ) + 2p: 1 cf{p s ,p7 1 )ci(p^p; 1 )-p; 1 [^(p i ,p7 1 )] + p7 1 cf (cf (p^pr^cf (p^pr 1 )) 
-c?(p-\ci( Pi ,p^))-ci(p;\c?{p t ,p^))+cf(p 7 \ [cf (P..P,- 1 )] 2 ) 

J •'J 3 a 3 tF J J 1 - ° 3 - 1 3 *- o 3 -i o £ 

+pr 1 [cf(p j ,pr 1 )] -cf^ScfCp.^-^-cffr-ScfCp,,?- 1 )) 

+2cf { P -\ [cfip^p; 1 )] 2 ) -cf (p r \cf (p,,^- 1 )) 

+2cf (p.-Scffc.p^cffc.p;- 1 )) -cf (p; 1 , [cf^.p; 1 )] 3 ) +cf (p 3 - 1 ,cf(cf(p 3 ,p 3 - 1 ),cf(p 3 ,p 3 - 1 )) 

= -p j - 1 cf(p 3 ,p; 1 )+p j - 1 [cf(p j ,pr 1 )] 2 -cf(pr 1 ,cf(p i ,pr 1 ))-cf(p j - 1 ! cf(p j! p j - 1 )). 

4.3 The functional calculus with the magnetic Schrodinger operator 

Using the results recalled in Section[3]we conclude that for any function g e Cq°(M), we can define by the functional 
calculus for self-adjoint operators a bounded operator g(H) that can be computed using Formula (12.31) . Then, using 

-— — B 

the fact that (H — ^) _1 = Dp h {r?) , we conclude that we can compute directly the symbol of g{H) =: Dp h (g(F) h ), 
9(F) l(X) = -7T- 1 km // ^(\,n)r^(X)d\d», VIe5, (4.4) 

defining a function of class S[ '~ (X). Using our previous Proposition 14.41 we obtain: 
Proposition 4.7 

,-^_- B 

The 'magnetic' symbol g{F) h admits for any K G N an asymptotic expansion in ft of the form 

gJf) B h (X)= J2 Vg?lF}(X) + h K g^{F}(h,X) 

0<j<K-l 

where the terms gf[F](X) only depend on the magnetic field B and its derivatives up to order j — 1 evaluated at 
X , and the rest g^ [F] (ft, X) only depends on the magnetic field B (in a non-local way). Moreover g B [F] € < S , 1 " J (X) 



tin 






dgB[F]£S?'->(X). 



Let us consider a function g s Cq°(R) such that its support E g is contained in ] — oo, ao[. As the spectrum of 
H is discrete in this region, and S g is a compact subset of ] — oo, Oo[, we deduce that g(H) is finite rank and thus 

~-^—B 

trace-class. Moreover, g(H) = Op (g{F) h ) is an integral operator having the integral kernel (see [ME]) 



K A [g(H)](x,y) := (e-*4,„l A ) #" 



9(F)n 



x + y 

—z—^-y 



with $ the inverse Fourier transform in the second variable (as defined on distributions on S = X x X'). In order 
to study the regularity properties of this integral kernel we shall use our formula (J4.4I) and write that 



where 



rfUPO = (F - (A + lf i))- 1 (X) + Xf+iJX) 



-n- 1 lim ff %(\^){F-{\ + iri)- l d\dn={goF)EC^(~) 

e^-0+ J J| M |> e oz 



and 

*?+* = rf +v - (F - (A + z/i))- 1 = rf +i/ jf (l - (F - (A + in))${F (A + z^))" 1 ) G Sf 2 - 2 (*) 
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due to our Proposition [33] applied to 1—(F— (A + i^))(jf (F — (A + i/i)) -1 . Let us remark that in 2 or 3 dimensions, 
£~ S^ 4 (X) is contained in the space of jointly continuous functions on X x X (by the Riemann-Lebesgue Lemma). 
Let us recall that for trace-class operators with continuous integral kernels we have the following property. 

Proposition 4.8 

Suppose T £ Bi(L 2 (A")) has an integral kernel K[T] £ C(X x X). Then the following limit exists and we have the 
equality 



lim / dxK[T](x,x) =TrT. 

i— >oo J 



\x\<R 



Let us discuss now the case d > 4. We come back to formula (|4.4[) and use Proposition 14.41 with the above 
observations. 

g~{F)l{X) = -it- 1 lim ff ||(A, M ) rf +iM (X)dAd/x 



, / 1(A,M) 
^+././t M |> £ dz 



e^0+ J J| M |> e OZ 



0<j<N-l 



(4.5) 



Taking now N > d and taking into account Proposition 14.41 we conclude that tjy(A + ifi; h) e S{ ' and by the 
Riemann-Lebesgue Lemma it has a continuous Fourier transform (with respect to the £ variable). Let us consider 
the main terms in the expansion (J43J) for N > d. Taking into account Proposition 14.51 we have to study integrals 
of the form 

lim / / gj(A^) f p {x^PllZ (*»fl - ff (*.0 (2 + p_ 1} , Jto || >e §(A^)(^ +1 ^)(^0 



o+././i„i> c 9 



= (2 ( + 1 p ) - +1 l)! fpB( "'° [( ^ +lg)0F]( "'°- 

We can evidently use again the Riemann-Lebesgue Lemma to obtain continuity of the Fourier transform (with 
respect to the £ variable). Putting all these results together we obtain the following statement 

Proposition 4.9 

For B a magnetic field with components of class BC°°(X) and H the self-adjoint operator defined in Proposition ^. 1\ 
if g G Cq°(R) has compact support S g C (— oo, do), then g{H) is a trace-class operatox\ and we have the formula 



Tig(H) = / dx3~ 
J x 

so that Tr g(H) only depends on the magnetic field. 



ml 



(./■.())= / dXg(F)l(x,0, 



4.4 End of the proof of the semiclassical trace formula 

4.4.1 Comparison of the theorems using Agmon estimates 

We are now ready to consider the semiclassical expansion of the trace formula starting from Proposition 14.91 and 
using the semiclassical expansions computed previously. Before doing that let us come back more in detail at the 
remark in HMR| that due to the exponential decay of the eigenfunctions (Agmon estimates |Ag| ) one can modify 
the potential outside a compact region by polynomially bounded terms with only an exponentially small change (of 
order exp{— c/h}) in the eigenvalues situated in any compact part of the discrete spectrum. A simple inspection of 
the proof in |Ag| shows that the same exponential decay estimates can be obtained for the magnetic Schrodinger 
operator so that we can apply the same arguments from [HMR to our 'magnetic' situation. Here is the basic 
proposition. 

Proposition 4.10 

Let (A, V) and (A, V) two pairs of electro-magnetic potentials satisfying Hypotheses \1.1\ and \1.2l Let E verify 

• E < min(Ey, £^>) , 

• U E :=V- 1 (]-oo,E[) = V- 1 Q-oo,E[), 

1 In fact it is even finite-rank. 
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• V — V on Ue and A = A on Ue , 

and let H and H the corresponding magnetic Schrodinger operators. Then for any g G Cq°(R), such that supp g C 
] — oo,E[, then Tr g(H) and Tr g(H) have the same semiclassical expansion modulo 0(h°°). 

It is enough to observe that, for any e > 0, the eigenfunctions corresponding to eigenvalues of H (resp. H) less 
than E — e ) decay exponentially in any compact outside of Ue—^- 

This can be used in the following way. 

Proposition 4.11 

Let (A, V) satisfy Huvotheses \l.l\ and \1.2\ and let E < T,y , then there exists a pair (A, V) such that the assumptions 
of Proposition \4-10\ are satisfied with in addition A and V bounded (with all the derivatives). 

The proof is easy. We can indeed consider a C°° increasing function \ on K such that 

x (t) =ton]-oc,-(E + £ y )[, X '(t) = on }-(E + 2S y ), +oo[. 

We can then take V = xOO- It is not difficult to modify A outside V _1 (] — oo, ^(E + Sy)[) to get a C°° bounded 
magnetic potential. 

As a consequence, it is enough for proving Theorem 11.31 to prove it with A and V of class C°° and bounded. 
Hence we can work at the intersection of the two calculi and use either the results of the Weyl's calculus or of the 
adapted magnetic calculus. 

4.4.2 The case with boundary 

Let us consider the case of the Dirichlet realization in a bounded open set $7, then it is easy to modify the comparison 
argument of the previous subsubsection in order to obtain the following theorem. 

Theorem 4.12 

Let A and V be C°° potentials on f2 and assume that 

miV{x) < inf V(x). 
ie!) x£dn 

Then, with H the Dirichlet realization of Pa in SI, there exists a sequence of distributions T B G 2?'(R), (j G N), 
such that, for any e > 0, for any N G N, there exists Cn and hjsr, such that if 

g G C£°(IR) , with suppg c] - oo, inf V - e[ , 

xEdfl 

then, : 

(2Trh) d Tr g(H) - ^ h j Tf(g) < C N h N+1 , \/h g]0, h N ] HI. (4.6) 

0<j<N 

More precisely there exists kj G N and universal polynomials Pg(u a , ua,j,fc) depending on a finite number of variables, 
indexed by a G N 2d and (3 G N d , such that the distributions: 

Tf(g)= E [9 (i) (F(x,0)Pe(dZ ti F(x,0,dPB jk (x))dxdt, (4.7) 

o<e<kj J 

Finally, T s = for j odd. 

Remark 4.13 The polynomials are the same as in Theorem \1.3\ . In particular they are independent ofQ. 

Using a (small extension of) the comparison proposition, one can modify the potentials in the neighborhood of 
dfl and then extend outside of Q without modify the asymptotic of Tr g(H) and then use the results obtained in 
the case of R d . 

Remark 4.14 Note that we have not done any assumptions on the topology of Q. Hence we have also that this 
expansion depends only on the magnetic field for cases where one can get various generating magnetic potentials 
which are not in the same cohomology class. 
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4.4.3 The odd coefficients vanish. 

To prove this result, one first observes that we have 



\h\ d Tr g(H) £ h j Tf(g) 

0<j<N 



< C N h N+l , \fh e [-h N , h N ] \ {0} . 



(the ft-pseudodiffcrential calculus can be extended to h < 0) and using the complex conjugation one obtains that 
the trace of g(H) is unchanged when h^r —h. Hence the odd coefficients are 0. 

References 

[Ag] S. Agmon: Lectures on exponential decay of solutions of second order elliptic equations. Mathematical notes 
of Princeton university n° 29 (f 982). 

[DiSj] M. Dimassi, J. Sjbstrand: Spectral asymptotics in the semi- classical limit. London Math. Soc. Lecture Notes 
Series 268, Cambridge University Press (1999). 

[Fels] B. Fedosov and T. Ispiryan: On perturbation of Morse Hamiltonians. Int. J. Geom. Methods Mod. Phys. 4, 
No 8, p. 1269-1283, (2007). 

[Hcl] B. Hclffcr: Introduction to the semi- classical analysis for the Schrodinger operator and applications. Springer 
Lecture Notes in Mathematics, n°1336. 

[Hc2] B. Hclffcr: h-pseudodifferential operators and applications: an introduction, Tutorial lectures in Minneapolis. 
The IMA Volumes in Mathematics and its applications, Vol. 95, Quasiclassical Methods, Springer Verlag (1997) 
p. 1-50. 

[HR1] B. Hclffer and D. Robert: Calcul fonctionnel par la transformer de Mellin et applications, Journal of 
Functional Analysis, Vol.53, No 3, oct. 1983. 

[HR2] B. Hclffer and D. Robert: Riesz means of bound states, semi-classical limit connected with a Lieb-Thirring 
conjecture, Asymptotic Analysis 3, p. 91-103, (1990). 

[HMR] B. Hclffer, A. Martinez and D. Robert: Ergodicite ct limite semi-classique, Commun. in Math. Physics 
109, p. 313-326, (1987). 

[HS1] B. Hclffcr, J. Sjostrand: Multiple wells in the semi-classical limit I, Comm. in PDE, 9(4), p. 337-408 (1984). 

[HS2] B. Hclffer and J. Sjostrand: Effct tunnel pour l'equation de Schrodinger avec champ magnetique, Ann. 
Scuola Norm. Sup. Pisa, Vol XIV, 4 (1987) p. 625-657. 

[HS3] B. Hclffer and J. Sjostrand: Equation de Schrodinger avec champ magnetique et equation de Harper, 
Partie I Champ magnetique fort, 

Partie II Champ magnetique faible, l'approximation de Peierls, 
Lecture notes in Physics, No 345 (editors A. Jensen and H. Holdcn), p. 118-198. 

[HS4] B. Hclffcr and J. Sjbstrand: On diamagnetism and the de Haas- Van Alphcn effect, Annales de 1'IHP, section 
Physique thcorique, (52), 1990, p. 303. 

[Hoi] L. Hbrmandcr: The spectral function of an elliptic operator, Acta Mathematica 121, p. 193-218 (1968). 

[Ho2] L. Hbrmander: Fourier integral operators I, Acta Mathematica 127, p. 79-183 (1971). 

[Ho3] L. Hbrmander: The Weyl calculus of pseudodifferential operators, Comm. Pure Appl. Math. 32, p. 359-443 
(1979). 

[Ho4] L. Hbrmandcr: The analysis of linear partial differential operators. Grundlehren der Mathematischen 
Wissenschaften, Springer Verlag (1984). 

[IMP1] V. Iftimic, M. Mantoiu, and R. Purice: Magnetic pseudodifferential operators, Publ. RIMS 43, No. 3, 
p. 585-624, (2007). 

[IMP2] V. Iftimic, M. Mantoiu, and R. Purice: Commutator criteria for magnetic pseudodifferential operators, to 
appear in Communications in Partial Differential Equations 2010. 

17 



[Iv] V. Ivrii: Microlocal analysis and precise spectral asymptotics. Springer Monographs in Mathematics (1998). 

[La] N. P. Landsman: Mathematical topics between classical and quantum mechanics, Springer- Verlag, New- York, 
1998. 

[LMR] M. Lein, M. Mantoiu, and S. Richard: Magnetic Pscudodifferential operators with coefficients in C*- 
algebras, preprint arXiv 0901.3704 and submitted. 

[MP] M. Mantoiu and R. Purice: The magnetic Weyl calculus. J. Math. Phys. 45 (4), p. 1394-1417, (2004). 

[MPR] M. Mantoiu, R. Purice and S. Richard: Spectral and propagation results for magnetic Schrodinger operators 
: a C*-algebraic framework. J. Funct. Anal. 250, p. 42-67, (2007). 

[Rob] D. Robert: Autour de I 'approximation semi-classique. Progress in Mathematics No 68, Birkhauscr (1987). 



18 



